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ABSTRACT 

In the present study, we analyze the effects of a flux of Alfven waves acting 
together with radiation pressure on grains as an acceleration mechanism of the 
wind of late-type stars. In the wind model we simulate the presence of grains 
through a strong damping of the waves, we use a non-isothermal profile for 
temperature coherent with grain formation theories. We examine the changes 
in the velocity profile of the wind and we show that if the grains are created in 
the region 1.1 < t/tq < 2.0 their presence will affect the mass loss and terminal 
velocity. The model is applied to a K5 supergiant star and for Betelgeuse 
(a Ori). 

Subject headings: MHD waves - dust - stars: late-type - mass loss 
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1. Introduction 



Mass loss from late-type stars is evidenced by the blueshifted circumstellar absorption 
lines present in their spectra. Giant and supergiant stars eject mass at a rate of 
~ 10^^'^ — 10~^ Mq yr^^ with terminal velocities Mqo ~ 10 — 50 km s~^, which is lower than 
the surface escape velocity {veo) (|Deutsch 1956| ; [Weymann 1962| ; [Habing 1996| ; [Lamers ^ 
ICassinelli 1999| ; [Carpenter et al. 1999|) . An outward-directed flux of Alfven waves has been 
proposed for driving these winds ( [Belcher fc Ulbert 1975| [Haisch, Linsky fc Basri 1980 



[Hartmann fc IVlacGregor 1980^ [Jatenco-Pereira fc Upher 1989a| [IVlacGregor fc (Jharbonneau 
1994 ; Uharbonneau fc Macgregor 1995| ). For late-type giant and supergiant stars, this 



nonthermal energy flux must be deposited in the region of subsonic flow to explain the 
observed u^o < feo ( [Hartmann fc MacGregor 1980[ ; [Leer fc Holzer 1980|) . 



We have direct evidence for Alfven waves in the solar wind. Waves are observed with 
large perturbations in the magnetic fleld, with negligible density perturbations, indicating 
the presence of Alfven waves. Since the early observation of MHD waves in the solar 
wind ( Pelcher fc Davis 1971|) , various authors have suggested that Alfven waves could 
be important in transfer momentum to the wind. We know from the solar coronal holes 
observations that the density is likely to fall off faster than (e.g. Hollweg 1981[ ; Poschek 
et al. 19"^ ), while the magnetic fleld strengths is proportional to ( Linsky 1987[ ). In 
our simulations we show that this behavior exists up to r < 3ro. The energy flux per unit 
area transported by the wave, 4>a, is 0a — e'^A — (l/2)po(<^'y^)'^A, where e is the energy 
density of the waves. This energy flux is constant, when there is no damping, or decreases 
due to damping, such that e decreases with r. Since the pressure associated with the wave 
is proportional to e, then this pressure also decreases with r. The result of this is a pressure 
gradient that accelerates the gas ( [Castor 1981 



Alfven waves have been used to explain the heating of stellar corona and the production 
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of stellar winds in many regions of the Hertzsprung-Russell (HR) diagram. They are used 
to explain the solar wind acceleration (|Jatenco-Pereira fc Opher 1989b| ; [Jatenco-Pereira 



Opher fc Yamamoto 199^ , cool supergiant stars (|Charbonneau fc Macgregor 1995 



Airapetian et al. 2000| ) and winds in hot stars (|Cassinelli 1979t [Underhill 1983| ; |dos 



Bantos, Jatenco-Pereira, fc Opher 1993|) . These waves can also be important in T Tauri 



variables ( [Hartmann, Edwards fc Avrett 195^ ; [Jatenco-Pereira fc Upher 1989c| ; [Vasconcelos 



Jatenco-Pereira fc Upher 200(J| ). The origin of Alfven waves was suggested by [Upher & 



Pereira (1986)[ to be the annihilation of twisted magnetic fields near the surface. These 



annihilation sites that could also explain the variable hot-spots observed in late-type stars 
(Gray 2000, 2001; [Young, Baldwin, Boysen et al. 2000|) . The waves can also be generated 



by convective motions inner in the envelope, that create perturbations in the magnetic 
fields. 

On the other hand, in cool supergiant stars, the presence of grains are well know 
determined by infrared excesses ( [Woolf fc Ney 1969| ). However, the region of nucleation 
is not precisely determined. These grains are supposed to be the main responsible for 
radiation absorption, and also are supposed to contribute in the acceleration of the wind 
(e.g. Jura 198^ ; Elitzur fc Ivezic 2001 ). The temperature profile of the wind for these 



stars is also not very well known, but by nucleation theories ( [Gail fc Sedlmayr 1987] ) the 



temperature range for grain formation is 400 K <T < 2500 K. 

Besides, the high luminosity of these stars have lead many authors to include the 
radiation pressure as a mechanism to explain the mass loss. MacGregor and Stencel (1992) 
studied the interaction between gas and dust in the curcumstellar envelope of late-type 
stars. In their analysis the grains and gas are well coupled if the grain radius is ~ 10^^ cm 
while for small grains < 5 x 10~^ cm, the collision momentum transfer is insufficient to 
induce expansion of the atmosphere whole. 
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Jatenco-Pereira and Opher (1989a) (hereafter JPO), presented a model for mass loss 
in a cool supergiant K5 star with M = IGMq, where a flux of Alfven waves flux {4>a) is 
the main mechanism for accelerating the wind. The model solves the equations along the 
magnetic field axis, and assumes (i) a divergent geometry for magnetic field lines, (ii) a 
pure plasma (without dust) and (iii) isothermal wind. Holzer, Fla and Leer (1983) argued 
that wind driven by Alfven waves from late-type stars only would be possible in an average 
ion-friction damping length ~ 0.85 — l.Oro. Since Lf is proportional to (where P 
is the average Alfven period), this requires that nature fine-tune P ~ 1.77 x 10"^ s. For 
long wave period (P > 10^ s), the wave damping length is large and the ion-neutral friction 
damping is not important (Pe Pontieu, Martens fc Hudson 2001| ). On the other hand. 



Davila (1985) emphasized that on the surface of the Sun one might expect that boundary 
effects may be important, and a WKB approximation may not be valid whenever the MHD 
wavelength (A) becomes greater than d ~ 10^° cm ~ 0.14ro, the characteristic transverse 
dimension of a solar coronal hole, in despite of Usmanov et al. (2000) that used a WKB 
approximation for the solar wind, in comparison with Ulysses data. We assume an Alfven 
wave flux with an averaged wavelength A < 0.5ro such that the WKB approximation is 
valid. JPO used a wave period P > 10^ s and discussed other possible damping mechanism 
for Alfven waves such as nonlinear damping, resonant surface damping and turbulent 
damping. They show that these damping processes can produce acceptable winds in 
late-type giant stars. We extend the JPO model including a simulation of a strong damping 
of Alfven waves, possibly due to a grain interaction, and the effect of radiation pressure in 
the wind acceleration. The changes in the velocity proflle of the wind are compared with 
the pure plasma and isothermal model. 

In section 2, we briefly describe the JPO model. In section 3, we present the inclusion 
of grain interaction with Alfven waves in a non-isothermal temperature proflle and the 
radiation pressure on grains. In section 4, we apply this magnetic-radiation model in a 
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K5 supergiant star and aOri and discuss the results. In section 5 we summarize our 
conclusions. 



2. The Model 

The model suggested by JPO for mass loss in late-type stars takes into account our 
knowledge of coronal holes in the Sun. From surface, the magnetic field configuration is 
divergent, and after a transition radius, tt, it becomes radial. The net non-radial expansion 
at this region is given by the factor of F = ^ = 7.26 (r < 3Rq) (|Munro fc Jackson 1977| ). 



JPO used F = 10 for supergiant stars. They used the diverging geometry for the magnetic 
field suggested by |Kuin fc Hearn (1982)| and Parker (1963)| given by 



A{r) = A{ro) 

Vro 

where A{r) is the cross section of the geometry at a radial distance r, tq is the initial radius 
and S* is a parameter that determines the divergence of the geometry up to the transition 
radius r^. The geometry is radial for r > tt- JPO analyzed the effect of the change in 
the opening angle of the magnetic field on mass loss and terminal velocity of a "typical" 
K5 supergiant star. They used a flux of Alfven waves as the acceleration mechanism of 
the wind, with non-linear, surface Alfven wave absorption and turbulent damping. In the 
model, for a given initial damping length at tq an upper limit of 5* exist when mass loss 
stops, Smax- They found acceptable solutions for 3 < 5 < 7. They showed for a range of 
initial damping lengths, Lq ~ 0.1 — 0.2ro, and an isothermal atmosphere, T ~ 10^ K, that 
a divergent geometry can produce the observed large mass loss rates and the small ratios 
Uoo/veo ~ 1/2 of the supergiants. 

The momentum equation of the wind may be deduced from continuity equation, 
conservation of momentum and action. The equation of motion (e.g., eq. [52] of |Holzer, Fla 
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& Leer 1983; eq. [20] of JPO) is 

Idu _ Z N 
u dr r D ^ 

where u is the ejected velocity at the radius r, Z = S for < r < and Z = 2 for r > r^. 
and D are given by 



D = u'- v' - - 



1 /I + 3Ma\ 2 



*^ 4 V 1 + 



j {^v') , (4) 



where = B / ^Anp, is the Alfven velocity, Vth = ykT/{jj,mH), is the thermal velocity, 
M4 = u/va, is the Mach number, u is the wind velocity, L is the damping length of wave 
and ^) the velocity associated to the magnetic field variation. 

The equations were solved in 1-D, along the magnetic field lines, in a cylindrically 
symmetric geometry, assuming the wave flux consistent with 6B/B << 1 and P > 10^ s, 
which is implicit in the initial damping length. In the present study the fluctuations in 
velocity, magnetic field and density are small allowing a linear treatment of Alfven waves 
( [Boynton fc 'Ibrkelsson 199(j| ), in despite of e.g. Ofman & Davila (1998) model for solar 
coronal holes, where the fluctuations of density and velocity were large and a nonlinear 
treatment of the Alfven waves were necessary. 



3. Grains simulation, temperature profile and radiation pressure 

3.1. Simulation of strong damping of Alfven waves 

In JPO model the Alfven waves were supposed to be damped by nonlinear process with 
damping length, Li, resonant surface damping, L2, and turbulent damping, L3, separately 
although they may contribute jointly. The damping length for each case is given by: 



- 8 - 



^10 



{dv^} \vaoJ 



(5) 



rJ \vaoJ 



Ma 



(6) 



f U + Va 



^30 



s 

r \ 2^ 



(7) 



Vuo + vaoJ \roJ \ {6v^) 

where Lio, L20, L-^q are the initial damping length for each mechanism and Uq is the initial 
plasma velocity. 

The above damping lengths were used in pure hydrogen plasma, i.e. there is no grain 
formation. Experimental ( [Elflmov et al. 1997 ), and theoretical ( Havnes, Hartquist fc 



Phillip 1989|) , researches show that Alfven waves are suddenly damped in "dust plasma" 
environments. Impurities in the wind may increase the dissipation of the Alfven waves. In 
the case of a stellar wind, inhomogeneous turbulent regions can occur immediately after the 
sonic point and they may manifest themselves as a sudden decrease in the damping length. 
We suppose that a sudden decrease in the damping length is correlated to the presence of 
grains although we do not treat them physically. This local dissipation may be simulated 
as an exponential decay of the damping length of the wave with distance given by: 



L grain = -^r ~ -^ri GXp , (8) 

where ri is the initial location in the wind where the grains can form which corresponds to 
the damping length L^^ and A is a constant that shows how suddenly is the damping decay. 

The determination of the grains formation region is highly uncertain, since this region 
can not be determined by the observations. Gail, Keller and Sedlmayr (1984), Gail and 
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Sedlmayr (1986, 1987) and Clayton et al. (1992, 1995) presented a theoretical nucleation 
model. Their results determines the region Itq < r < 2ro for SiO grains formation at cool 
supergiant envelopes, that is named here by "grain formation point" (GFP, hereafter). 

Analyzing the dependence of the three previous damping mechanisms in the JPO 
model vs. r (eqs. 5, 6 and 7) we note that all of them have the same behavior in the inner 
regions of the wind. They initially rise and peak for r < 2ro as we can see in Figure |l]. 
It means that considering the region of nucleation given by Gail, et al. (1984), the model 
is independent of which initial damping mechanism we use, since differences only occur 
at large distances. For comparison we plot in Figure 1 these three damping lengths as a 
function of r and the simulation of strong damping of Alfven waves given by eq. 8. In our 
simulations we used the nonlinear mechanism (Li) from the stellar surface until the region 
of grain formation and thereafter the wave damping assumes the form given in eq. 8. 

[EDITOR PLACE FIGURE 1 HERE] 



3.2. The temperature profile 

Observations of the solar atmosphere show that the temperature (~ 6000 K) decreases 
reaching a minimum temperature at about 500 km above surface and increases slightly 
after, to ~ 10'^ K, into the chromosphere, and then suddenly jumps through a narrow 
transition layer to ~ 2 x 10^ K in the corona. The solar wind is an extension and a 



super-sonic expansion of the solar corona into interplanetary space ( [Priest 1985| ) 



As in the Sun, we assume that the same behavior of temperature profile takes place in 
the giant and supergiant stars. This profile must be consistent with grain formation, since 
the grains can be formed only at low temperatures regions. 

In this work we introduce a temperature profile from a minimum ~ 2000 K, above the 
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photosphere, similar to the profile obtained by Hartmann & Avrett (1984) and a rapidly 
increase to ~ 10^ K ai r = 2ro. This profile is in agreement with the Gail and Sedlmayr 
(1986, 1987) theory. Rodgers and Glasshold (1991) modelled the thermal equilibrium of the 
alpha Ori atmosphere, assuming that heating occurs by grain-gas collisions and by radiation 
absorption and the cooling occurs by adiabatic expansion and radiative de-excitation. They 
found a decreasing profile for r > 2ro, from 10^ to 2000 K aX r = lOrg. By this way, we 
combined the solar observations, for sudden increase of temperature in chromosphere, to 
the Rodgers & Glasshold (1991) model. 



3.3. Radiation pressure 

Many authors have tried to explain the mass loss processes of luminous stars like 
Wolf-Rayet using only radiation pressure as the accelaration mechanism. The great problem 
that arises in this star is called "the momentum problem" , that consists of a radiation field 
momentum lower than the momentum of the wind ( Abbott 1978| ): t] = Muoo/(-^*/c) is in 



the range 3 — 30 ( [Barlow, Smith fc Willis 198 1|) . dos Santos, Jatenco-Pereira & Opher 
(1993) presented a model for mass loss in Wolf-Rayet stars where a fiux of Alfven waves 
acting together with radiation pressure has enough momentum fiux and energy fiux to drive 
the wind. They assumed a plasma without dust and isothermal wind. 

In order to verify if radiation pressure can be important in driving late-type star wind 
we include the radiation pressure into the wind equation using an effective escape velocity. 
Taking a pure radiative driven wind, the force exerted on a grain of radius a situated at a 
distance r from a star of luminosity and mass is: 
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where Q is the mean extinction efficiency factor of the grain, and a ~ 10 ^ cm is the grain 
radius value. 

Using the Newton's Law, the wind equation becomes: 



dr p dr p ' 



or 



du IdP GM^, , , , 

+ + ^ 1-r = 0, 11 
dr p dr r^ 

where the factor F is given by: F = ^ I "cGaT ' "^i^^ y the ratio of the grain number density 
to the gas mass density. 

If we caU: 



-^(-e.er = -^(l-r), (12) 

we define the effective escape velocity.. 

The escape velocity (12) is introduced in equation (3) in the region where the grains 
may affect the wind dynamics i.e. for r > GFP (with l.lro < GFP < 2ro). We assume that 
most of the radiation is absorbed by the grains in this region. This radiation absorption is 
partially responsible for grain destruction and acceleration. These accelerated grains are 
dynamically connected to gas ( |Lafon fc Berruyer 1991| ; [Liberatore, Lafon fc Berruyer 2001| ), 



so that the grains transfer momentum to the gas inducing the mass loss. 



4. Results 



This magnetic-radiation driven winds model consider a steady, cylindrically symmetric 
and 1-D wind emanating from a divergent magnetic field geometry with a temperature 
gradient similar to the Sun atmosphere. We assume that a flux of Alfven waves, which 
properties were defined in earlier sections, are generated near the stellar surface and 
propagate outward. Since in the inner region of the wind the damping mechanisms already 
studied in JPO have the same behavior we so use the non-linear one untill the GFP, after 
that we use the damping length given by eq. 8. In the region of grain simulation (r > GFP) 
the radiation pressure in grains, as well as the grain damping, are taken into account from 
GFP to 300ro. We applied the model to the K5 supergiant, (model 6 of Hartmann & 
MacGregor (1980)) and aOri. 



4.1. K5 supergiant star 

We used M, = IGM©, 5 = 10 G, tq = 4OOi?0, = 3 x lOU©, 
(pAo = 3.36 X 10^ erg cm~'^ F = 0.3, S* = 5, A = 0.1, the initial damping 
length Lq = 0.2ro and the density on wind basis p = lO^^^gcm^^. The typical velocity and 
mass loss rate observed for these stars are of order Uoo ~ |feo and M = 5 x 10^^ Mq yr^^ 
(see e.g. [Gharbonneau fc Macgregor 1995| ; [Whitelock et al. 1995| ; [Lamers fc Gassinelli 1999 



[Lindqvist et al. 2000|) . 

First we performed some simulations in an isothermal wind in order to evaluate the 
influence of the grains on the terminal velocities, found by previous pure plasma model. The 
introduction of damping by grains has an important effect in decreasing the flux of Alfven 
waves much more effectively than the others damping mechanisms. As a consequence the 
terminal velocity of the wind reaches lower values as compared with JPO model. Since 
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the grain formation region is not well known, we introduced the damping by grains at 
different distances from the stellar surface. The models show that as closer to the surface 
the grain if formed, lower is the asymptotic velocity. This effect is caused by the lose of 
wave energy by the strong damping, which leaves the wind unsupported against gravity. 
Decreasing the non-radial expansion factor 5* we can coUimate the wave flux into a smaller 
area accelerating the gas more efficiently. On the other hand, we found that for r > 2ro the 
presence of grains is not important for the wind dynamics. The properties of the model can 
be seen in Table 1. 

[EDITOR PLACE TABLE 1 HERE] 

The velocity curve obtained by solving the equation of the wind (2) is plotted in 
Figure |^. For comparison we plot the velocity profile for a pure plasma and isothermal 
wind (dotted line), and for the magnetic-radiation model (solid line). We may note that 
the terminal velocity for the magnetic-radiation model is higher than the JPO concluding 
then that the introduction of grain and temperature profile affect the wind dynamics. The 
grain presence diminishes the asymptotic velocity, while the introduction of a temperature 
gradient increases the velocity. For a given density, temperature, magnetic field and GFP 
(which are all observable) we may adjust the model to the expected values of u/veo and 
M, by varying the Alfven wave flux. The best fit occurs for an initial Alfven wave flux of 
5 X 10® erg cm~'^ s~^. On the figure 2, for a fixed GFP, we note that the terminal velocity 
is higher for lower S values due to the consideration of a radiative and non-isothermal 
atmosphere. In this sense increasing the S value we should have a decrease in the terminal 
velocity. However, for higher S values there was no sonic point. The observations can be 
reproduced if we increase the initial Alfven flux to 5 x 10® erg cm~^ s'^. In the flgure we 
mark the region where we assume the formation of grains. 

[EDITOR PLACE FIGURE 2 HERE] 
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Although for supergiant stars the momentum problem is not present, it is interesting to 
estimate the value of rj. For example, for K5 supergiant star using the above values we have 
7] ~ 0.4, lower than 1, but still high considering normal absorption coefficients. Using the 
above values for the K5 supergiant star, we can estimate the ratio {^) of the initial energy 
flux of Alfven waves to the radiation energy flux: ^ = (^aq/ {L^/ {ATtr^)) ~ 10^^. This value 
indicates that the energy flux of Alfven waves is only a few percent of the stellar luminosity. 
This relation may be used to estimate the Alfven wave flux on other stars like Betelgeuse, 
for example. 



4.2. Betelgeuse {aOri) 

Betelgeuse observations give the terminal velocity, Mqo ~ 15 km and a mass loss 
rate M ~ 2 - 5 x 10"^ Mq yr'^ (|Loup et al. 1993| ; |Knapp fc Woodhams 199l|) . We 



applied the magnetic-radiation model to this star with a temperature proflle given by 
Rodgers and Glasshold (1991). We used the following initial conditions: M^, = 8Mq, 
B = 10 G, ro = 3OOi?0, = IO^Lq and T = 0.7, S = 5, Lq = 0.2ro and the GFP 
located at r = 1.24ro. The Alfven wave flux can be determined by ^ = 10"'^. We found 
0AO = 6.9 X 10^ erg cm~'^ s^^. As in the previous section, we performed the calculations 
until r = 300ro. However, we show in Figure 3 the resultant velocity proflle until r = IOtq 
since for r > lOro the curve becomes flattened showing the same terminal velocity. 

[EDITOR PLACE FIGURE 3 HERE] 

We obtained Uoo ~ 16 km and M = 4 x 10^^ Mq yr~^, very similar to the 
observations. Betelgeuse also has high UV emission at r ~ 2ro region. This emission may 



be connected to the MHD waves damping ( pilliand fc Dupree 1996|) and with a sudden 



heating of this region, which is in agreement with our model. 
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5. Conclusions 

In this work we present a model to study mass loss processes in late-type stars, using an 
Alfven waves flux added to radiation pressure as accelerating mechanisms. We simulate the 
presence of grains in the sense of a strong damping of Alfven waves in the region r > GFP 
(which is about l.lro < GFP < 2ro). In this region we added the radiation pressure as an 
acceleration mechanism of the wind. In accordance with grain formation theories and in 
agreement with observations, we took into account a temperature profile for the wind going 
from ~ 2000 K above the photosphere, with a rapidly grown to ~ IQl^K at r = 2ro and a 
decreasing profile for r > 2ro. The models show that as closer as to the surface the grain 
is formed, lower is the asymptotic velocity and that the wind velocity is not changed if the 
region of grain formation is located above 2ro. Although a number of hypotheses have been 
made this hybrid model can lead to a more realistic wind and gives more constraint in 
the physical parameters of the wind. The model was applied to a K5 supergiant star and 
Betelgeuse, resulting in values that agrees with observations. Further studies will include 
the physical parameters of the grains as well as their distribution in the wind. 
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partial support, respectively. D.F.G. would like to thank the Brazilian agency FAPESP for 
providing computational equipment under grant (No. 96/00677-3). The authors would like 
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Figure captions 

Fig. 1. — The behavior of damping length as a function of the distance for resonant damping 
(sohd curve), nonhnear damping (dashed curve), turbulent damping (dashed-dot curve) and 
by grains (dotted curve). The initial damping length is Lq — 0.2ro for the three cases. For 
damping by grains (eq. 8) we used ri = 1.32ro, Lri = 0.75ro and A = 0.1. 

Fig. 2. — Velocity profile for JPO homogeneous and isothermal model (dotted hne) and 
for the hybrid model including grain presence, non-isothermal atmosphere and radiation 
pressure (solid line). The grain formation region is marked in the figure. 

Fig. 3. — Velocity profile as a function of distance for aOri. In this non-isothermal model 
we used a flux of Alfven waves added to radiation pressure as an acceleration mechanism for 
the wind. 
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